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Abstract 

In hadronic decays of Z bosons recorded with the OPAL detector at LEP, events containing 
b quarks were selected using the long lifetime of b flavoured hadrons. Comparing the 3-jet 
rate in b events with that in d,u,s and c quark events, a significant difference was observed. 
Using C?(a§) calculations for massive quarks, this difference was used to determine the 
b quark mass in the MS renormalisation scheme at the scale of the Z boson mass. By 
combining the results from seven different jet finders the running b quark mass was 
determined to be 

m h (m z ) = (2.67 ±0.03 (stat.)^ (syst.) ±0.19 (theo.)) GeV . 

Evolving this value to the b quark mass scale itself yields m b (m b ) = (3.95±o;(l) GeV, 
consistent with results obtained at the b quark production threshold. This determination 
confirms the QCD expectation of a scale dependent quark mass. A constant mass is ruled 
out by 3.9 standard deviations. 
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1 Introduction 



In Quantum Chromodynamics (QCD) the renormalisation group equation (RGE) governs 
the energy dependence of both the renormalised coupling as and the renormalised quark 
mass m q . The RGE for an observable R calculated for massive quarks q and measured 
at a scale Q, states that R is independent of the renormalisation scale \i 
expressed by 



which is 



_d_ 



n, \ 9 , . d 

+ p{a s )- l{a s )m q - — 

das om c 



R(Q 2 / f i 2 ,a s ,m q /Q) = 



with the (3 function /3(«s) of QCD and the mass anomalous dimension 7(03)- This 
equation can be solved by introducing both a running coupling constant as(Q 2 ) and a 
running quark mass m q (Q 2 ). In particular, the scale dependence of the b quark mass in 
the MS renormalisation scheme, mb(Q 2 ), is to four-loop accuracy given by 



m h (Q 2 ) 



'ois{Q 



2 x \ 12/23 



7T 



1 + 1 . 175 . (!*m) + L500- (!*m' + 



0.172 




7T 



2\4 



(2) 



taking the renormalisation group invariant^] mass rh-b as a reference, see e.g. 0. Analogous 
to a$(Q 2 ), an absolute value for m h (Q 2 ) is not predicted by QCD. A b quark mass of 4.2 
GeV measured at the production threshold corresponds to a running mass m h (Q 2 ) of 
about 3 GeV in interactions at the scale of the Z mass. The experimental observation of 
this running of the quark mass constitutes an important test of QCD. 

Studies of the flavour dependence of the strong coupling constant observed a difference 
in jet rates and event shapes between b events and light quark events, see e.g. This 
apparent deviation of a few percent from a flavour-independent coupling constant can 
be explained by effects of the large b quark mass. Second order matrix elements that 
have been calculated recently, taking finite quark masses fully into account [@-[7] can 
explain these experimental observations. Flavour independence of the strong interaction 
is a fundamental property of QCD. Assuming it holds, the second order matrix elements 
for massive quarks can be used to determine the b quark mass at scales different from 
production threshold. 

For the determination of the running b quark mass, the ratio of 3-jet rates in b events, 
-R3, over 3-jet rates in light quark events, i?3 USC , 



R du SC 



(3) 



1 i.e. independent of [i 
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has been proposed in 0]. This ratio is sensitive to mass dependent differences in gluon 
radiation from b and from light quarks. This or equivalent methods have been used in 
determinations of the b quark mass by DELPHI ||, ALEPH M and Brandenburg et al. 
using SLD data 0. 

In this paper we present a determination of the running b quark mass based on the 
variable B% using the large statistics sample collected with the OPAL detector at the e + e~ 
collider LEP at centre-of-mass energies close to the Z mass. Events containing b hadrons 
were tagged by identifying their displaced decay vertices. The light and the b quark 
contributions were deduced from the tagged and inclusive samples by a simple unfolding 
technique which relies only on the tagging efficiencies and fake tagging rates which were 
estimated by studying Monte Carlo events. 



2 The OPAL detector, data, and Monte Carlo simu- 
lation 



A detailed description of the OPAL detector can be found elsewhere |TT| . For the present 
analysis only e + e~ collisions collected in 1994 were included, as these provide sufficient 
statistics with a uniform detector configuration. The silicon strip micro-vertex detector, 
all central tracking detectors and the electro-magnetic calorimeter were required to be 
fully operational. Standard criteria for high multiplicity hadronic events |12|,[13|, which 
rely on a minimum number of measured tracks in the central tracking system and clusters 
in the electro-magnetic lead glass calorimeter were applied. The remaining background, 
mostly from two-photon processes and r-pair events, was estimated to be 0.5 % and 0.11 %, 
respectively 
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The silicon micro-vertex detector [|14[] is used for the identification of b quark events 
based on lifetime information. To account for the limited polar angle acceptance of the 
silicon micro-vertex detector operating in 1994, only events whose thrust vector pointed 
to the central part of the detector, | cos thrust | < 0.75, were considered^. After this cut 
about 924000 events remained. These events defined the inclusive sample. Tracks of 
charged particles recorded in the tracking detectors and clusters of energy recorded in 
the calorimeters were used for jet finding and were required to satisfy a set of standard 
quality cuts which are detailed in JT5| ]. The energy of each cluster associated to a track 
was corrected for double counting of energy using the momentum of that track [16| . 



To determine the efficiency and purity of the event selection and to correct for distor- 
tions due to the finite acceptance and resolution of the detector, about 4 million hadronic 



decays of the Z were generated by the JETSET program version 7.4 |17|], tuned to describe 
OPAL data |TB|. The generated events were passed through a detailed simulation of the 



OPAL detector [19] and reconstructed using the same procedures as for the data. The b 



2 OPAL uses a right handed coordinate system with the z axis pointing along the electron beam 
direction and x towards the centre of the LEP ring. The polar angle 9 is measured with respect to the z 
axis. 
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quark events in the Monte Carlo sample were reweighted to correspond to the most recent 
estimates for the parameter of the Peterson et al. fragmentation function 



the mean charged particle decay multiplicity in b events, ^charged EH' an< ^ ^he mean b 
hadron lifetime, tb |[22|| . The c quark events were reweighted to correspond to a recent 
estimate of the Peterson fragmentation parameter, e c |[L8| . The reweighting procedures 
are described in p3 
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3 Selection of b quark events 



The silicon micro- vertex detector was used in addition to the tracking detectors to measure 
the decay length of b flavoured hadrons in hadronic Z decays. The decay length is defined 
by the distance between the reconstructed primary vertex and the identified b hadron 
decay vertex. These vertices were reconstructed using the algorithm described in EM . In 



the procedure a cone jet algorithm |25j is applied to search for jets in each event, using 
a cone half angle with R = 0.55 rad and a minimum jet energy of 5.0 GeV, as in ||23|| . 
A common secondary vertex was searched for in such a jet by iteratively excluding the 
track with the largest x 2 contribution and repeating the fit until all x 2 contributions were 
smaller than 4. A minimum number of three tracks was required to form a vertex. For 
each event the vertex with the largest decay length significance L/<jl was determined, 
where L is the decay length and <jl its uncertainty. 

The distribution of L/<jl is shown in Figure |l|. A good agreement between data and 
simulation is observed in the region Ljo^ > 5, used to select a sample enriched in b events. 
In the Monte Carlo b quark events are selected with an efficiency of e = (70.23 ± 0.02)%, 
whereas the fake tag rate from light quark events being mis-identified as a b candidate 
is / = (7.20 ± 0.01)%, where the uncertainties are statistical only. 21.3% of all hadronic 
events are tagged as b candidates in the data, compared with 21.0% in the simulation. 
This small deviation will be discussed in section EO. 



4 Measurement of B- 



For the determination of -B3 jets were reconstructed using the standard JADE algorithm, 
its variants E, E0, P and P0 and the DURHAM and the GENEVA algorithms, all de- 
scribed in [ 2~6| and [ 2~7 1 . In addition, the CAMBRIDGE algorithm was used. However, 



since there is not yet a second order calculation compatible with our definition of -B3 for 
this jet finder available, it was not used to determine the b quark mass. All these jet find- 
ers combine the two objects with the smallest distance as measured using the distance 
measures in Table [I]. These two objects are combined according to the prescription given 
in the third column of Table [I] to form a new object. This procedure is repeated until all 
distances between objects are larger than a resolution parameter y cut . The number of jets 
is then given by the number of remaining objects. To limit the uncertainty related to the 
choice of a specific jet finder, seven different jet finders were used in this analysis. 
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Figure 1: Distribution of the decay length significance L/a L for data (points) and 
simulation (histograms). The contributions from the d, u or s and c quark events in 
the simulation are shown by different hatching, b quark events are shown by the open 
histogram. The vertical line indicates the cut chosen to select b quark events. The 
statistical uncertainty is also shown. 
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DURHAM 
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ft = Pi + ft 


GENEVA 


8£*£y(l - cos%)/9(£i + Ej) 2 


Pk = Pi + pj 


CAMBRIDGE 


2(1 - cos%) 
soft freezing if y {j > y cut 
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Table 1: Definitions of the distance measures and recombination prescriptions for the 
different jet finders used in this analysis |26|,|27|. ft, p i; Ei describe the four-, the three- 
momenta and the energy of particle i. 8ij is the angle between the three-momenta of 
particles i and j. The recombination prescription of the CAMBRIDGE jet finder is 
described in more detail in IEB . 



The double ratio B 3 (y cnt ) was determined from the 3-jet rates in the event sample 
enriched in b quarks and in the inclusive event sample at values for y cnt at which the 
predictions were calculated, see Table As the contribution of events not originating 
from hadronic Z decays is very small, the inclusive sample can be decomposed into a b 
quark sample plus a light quark sample. The number of events, A, in the inclusive sample 
and the number of events, T, in the b enriched sample can be written in terms of the 
number of b quarks events, A^ b , and of light quark events, N dusc : 

A = N h + N dusc (4) 
T = eN h + fN dusc (5) 

where e is the Monte Carlo tagging efficiency for the b quark events and / is the fake tag 
rate for light quark events. These two equations can be solved for A^ b and A^ dusc . A similar 
decomposition is valid for the number of 3-jet events, a 3 (y cut ), in the inclusive sample and 
the number of 3-jet events, t 3 (y cut ), m the b enriched sample, with the number of 3-jet b 
quark events, n^(y cut ), and the number of 3-jet light quark events, n dusc (y cut ): 

a 3 (y C nt) = n h 3 (y cut ) + n dusc (y cut ) (6) 
^(z/cut) = e 3 (y cut )nl(y cut ) + f 3 (y cut )n dusc (y cut ) (7) 

The variables e 3 and f 3 , which depend on y cu t, are the corresponding efficiency and fake 
tagging rate for 3-jet events. Solving these last two equations for n 3 (y cnt ) and n dnsc (y cu t) 
and using the similar equations for A^ b and A^ dusc leads to the following relation: 



B 3 (y C ut) — C h ad(z/cut) • C dct (y cut ) ■ 

= Chad(?/cut) • Cdet(?/cut) 
= Chad(?/cut) • Cdet(|/cut) 



R h 3 (y c 



Rt sc (y C ut) 

hi 



n h 3 (y cut )/N h 
n dnsc (y cn t)/N^ 
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Vcut 


^Dct 


^Had 


JADE 


0.02 


0.965 


1.016 


DURHAM 


0.01 


0.961 


0.989 


JADE EO 


0.02 


0.993 


1.073 


JADE P 


0.02 


0.978 


1.015 


JADE PO 


0.015 


0.985 


1.022 


JADE E 


0.04 


0.997 


1.049 


GENEVA 


0.08 


0.971 


1.031 


CAMBRIDGE 


0.01 


0.971 


1.017 



Table 2: Correction factors for each jet finder for detector distortions and hadronisa- 
tion effects. The JETSET Monte Carlo program was used to estimate the size of the 
corrections. 

hiVcut) ~ .Ml/cut) • «3(?/cut) T - e ■ A 

T-f-A t 3 (y cut ) - e 3 (y cnt ) ■ a 3 (y cut ) 

Since B 3 is a ratio, common correction factors for the individual 3-jet rates of b quark and 
light quark events cancel. We apply bin-by-bin correction factors for detector distortions, 
Cdet(?/cut), and hadronisation effects, Chad(2/cut), as shown in Eq. (||). The correction 
factors are defined as the ratio of the double ratio of the 3-jet rates for b over light quark 
events determined from the simulation at either the hadron or parton level divided by 
the same double ratio at detector or hadron level. The hadron level consists of particles 
generated by the Monte Carlo program with a mean lifetime greater than 300 ps. The 
partons which are present at the end of the parton shower in the generator define the 
parton level. The parton shower in JETSET, which we used to estimate the size of the 
hadronisation corrections, terminates when partons reach virtualities below a cut-off Qo, 
set to 1.9 GeV in the standard analysis. 

Figure |2| shows the y cut dependence of the correction factors for the DURHAM and 
the JADE E0 algorithms. The DURHAM scheme has the smallest, the JADE E0 scheme 
the largest hadronisation corrections of all schemes used in this analysis. The correction 
factors for the other jet finders are summarised in Table 

The efficiency and the fake tag rate for 3-jet events in b or light quark events, respec- 
tively, shown for the DURHAM and JADE E0 jet finder in Figure |3], depend slightly on 
the chosen y cnt value because of the difference in kinematics induced by the presence of a 
highly energetic gluon emitted at large angle. 

In Figure |] the measured B 3 ratio is shown both before applying corrections and after 
being corrected to the parton level. The detector correction factors, which take into ac- 
count kinematic biases induced by the b tagging, are usually larger than the hadronisation 
correction, which includes known decays of b flavoured hadrons in the b quark sample. 
Only for JADE E0 and JADE E these corrections are larger than the detector corrections. 
All other hadronisation correction factors are smaller than about 3%. 
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Figure 2: In a) and b) the correction factors of the B 3 ratio for detector effects, Cd e t, 
and in c) and d) for hadronisation effects, Chad, are shown as a function of y cut for the 
DURHAM (left column) and JADE E0 (right column) jet finders. The error bars show 
the statistical uncertainty of the Monte Carlo. The dashed vertical lines indicate the y cut 
values chosen to determine the b quark mass. 
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Figure 3: a) and b) show the efficiency e 3 (y cut ) for b quark 3-jet events for the DURHAM 
(left) and JADE E0 (right column) jet finders. Overlaid as a dotted line is e, the tagging 
efficiency for all b quark events; c) and d) show the fake tag rates /^(l/cut) fo r light quark 
3-jet events for both jet finders. Overlaid is the fake tag rate / for all light quark events. 
The error bars describe the statistical uncertainty of the Monte Carlo. The dashed vertical 
lines indicate the y cut values chosen to determine the b quark mass. 
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Figure 4: The uncorrected a) and b), and the corrected c) and d) B 3 ratios for the 
DURHAM and the JADE E0 jet finders. The uncertainties include the effect of finite 
statistics in both data and Monte Carlo. There is a substantial correlation between adja- 
cent bins. The dashed vertical lines indicate the y cut values chosen for the determination of 
the b quark mass. The bands in c) and d) display the total systematic and the statistical 
uncertainty added in quadrature. 
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5 Determination of 



To determine the b quark mass from the measured B3 ratio, parametrisations of the 
QCD predictions for the 3-jet rate for b and light quarks were used. For b quarks the 
predictions 0,|3O| cover masses in the range from 2 to 4 GeV in steps of 0.2 GeV at one 
specific y cut value for each jet finder. These y cut values are listed in Table [2|. For light 



quarks the parametrisation of [571] was adopted except for the DURHAM jet finder, for 



which a more recent parametrisation from p8f was used. 



Using these parametrisations the double ratio B 3 = R^/Rf 130 was derived: 



B 3 (y C ut,rn h ) 



# 3 dus %cut) 
(g) A b (y cutj m h ) + (f ) 2 B h (y cuu m h ) 

(f ) A(y mt ) + (f ) 2 [B(y cut ) - 2 ■ A(y cnt )} 



(9) 



where the coefficients A h (y cnt , rrib) and B h {i/ cnt) m\>) and A(y cut ) and B(y cnt ) parametrise 
the 3-jet rates for massive and massless quarks. Note that the B coefficients for massive 
and for massless quarks differ in their definitions because the 3-jet rate for massive quarks 
is normalised to the total bb(g) cross section, whereas the 3-jet rate for massless quarks 
is normalised to the hadronic born, i.e. only qq, cross section. To obtain a prediction in 
a finite order of as, the denominator of the double ratio was expanded in a Taylor series 
after cancelling one order of as in numerator and denominator. Due to this cancellation 
of one order of as, the QCD predictions in this expanded expression for B3 are of 0(as). 
Due to the additional dependence of B3 on the b quark mass, the dependence on the 
renormalisation scale x M = fi/Q enters in first order of as, as can be seen from the 
expanded expression: 

B 3 ( Vmt ,m h ,x») = A{y ^ +-• 

( B b (y cut ,m h ) A b (y cut ,m h ) B(y cut ) - 2 • A(y cat ) 



A(y cu t) A(y cut ) A(y, 



cut J 



2vT7o^b T7 r log a; L 10 

A{y C ut) 

where 70 is the first coefficient in the perturbative expansion of the anomalous mass 
dimension. The renormalisation scale parameter was set to unity at the renormalisation 



scale of 91.2 GeV and as was set to its world average of 0.1184 [31]. To parametrise the 



mass dependence of Eq.([TO|) in the range from 2 to 4 GeV a parabolic function 

Bl(m h ) = a l + 4 • ml , i = JADE, E, E0, P, P0, DURHAM, GENEVA (11) 

was used, where the coefficients and their uncertainty due to a finite Monte Carlo inte- 
gration sample for R% are given in Table ^|. The last column gives the x 2 /d.o.f. of the fits. 
In [[l(J a different parametrisation was chosen. The difference between this parametri- 
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2/cut 






a\ x 1000 


Aa\ x 1000 


rVd.o.f. 


JADE 


0.02 


0.9909 


0.0013 


-2.505 


0.125 


0.233 


DURHAM 


0.01 


0.9918 


0.0025 


-3.771 


0.230 


0.383 


JADE EO 


0.02 


1.0142 


0.0013 


6.093 


0.124 


3.283 


JADE P 


0.02 


0.9824 


0.0013 


3.845 


0.126 


1.580 


JADE PO 


0.015 


0.9756 


0.0021 


5.749 


0.186 


0.868 


JADE E 


0.04 


1.0043 


0.0021 


9.886 


0.211 


0.979 


GENEVA 


0.08 


1.0171 


0.0013 


-2.406 


0.132 


0.364 



Table 3: Coefficients for Eq. ( |XT| ) for the different jet finders. 



sation and the one we used in the relevant mass region of 2 to 4 GeV is smaller than 
the statistical uncertainty of the Monte Carlo integration. To calculate the derivative 
dA h (y cnt ,m\ } )/dmb in Eq. (|10"D the same parametrisation as in Eq. (|ll|) was fitted to the 
A h (y cu t, Trib) coefficient and the derivative was determined. 

In Table |] the measured B 3 values are shown, together with the results for obtained 
using Eq. ([TTJ) along with their uncertainties, described in more detail in the next section. 
Correlations between the different event samples entering the unfolding in Eq. (H) were 
taken into account. These correlations were calculated from the full Monte Carlo sample. 





2/cut 


B 3 


m b [GeV] 


JADE 


0.02 


0.9752 ±0.0048 ±0.0129 


2.51 ±0.42 ±0.99 ±0.16 


DURHAM 


0.01 


0.9532 ±0.0056 ±0.0166 


3.20 ±0.24 ±0.64 ±0.18 


JADE E0 


0.02 


1.0657 ±0.0035 ±0.0164 


2.91 ±0.10 ±0.50 ±0.26 


JADE P 


0.02 


0.9938 ±0.0040 ±0.0152 


1.73 ± 0.34 ± 1.93 ±0.11 


JADE P0 


0.015 


1.0140 ±0.0035 ±0.0129 


2.58 ±0.12 ±0.45 ±0.09 


JADE E 


0.04 


1.0649 ±0.0032 ±0.0154 


2.47 ±0.07 ±0.32 ±0.25 


GENEVA 


0.08 


0.9817 ±0.0094 ±0.0257 


3.84 ± 0.55 ± 1.27 ±0.12 


CAMBRIDGE 


0.01 


0.9717 ±0.0062 ±0.0192 





Table 4: For each jet finder the double ratio B 3 and the corresponding value for the b 
quark mass are given. For the double ratio B 3 the statistical and systematic uncertainties 
are given. For the b quark mass Mb the uncertainties listed are the statistical, the total 
systematic and the total theoretical uncertainties. 



6 Systematic and theoretical uncertainties 

Various sources of systematic uncertainty were investigated to assess their impact on the 
measured b quark mass. Selection criteria and parameter values in the Monte Carlo 
simulation were changed from their defaults and the entire analysis was repeated. The 
deviation of the mass value from the standard result was taken as a systematic uncertainty. 
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The investigations can be grouped into three classes according to the corrections and 
efficiencies in Eq. (§). These classes are either related to (i) detector and b tagging, (ii) 
the hadronisation correction uncertainty for the Monte Carlo generator, or (iii) theoretical 
uncertainty. For any systematic variation which has a positive and a negative contribution 
the larger of both was taken as the symmetric uncertainty. For the variations which have 
only one deviation, the varied cut on the decay length significance, the detector simulation 
and the modelling of the b quark fragmentation, the deviation was taken as the symmetric 
uncertainty. 

Tables |5] and |B| summarise the results of these checks along with their uncertainties 
assigned. The fairly large spread between the numerical values for the systematic varia- 
tions for the different jet finders is caused by different slopes for the parametrisation of 
the double ratio in terms of the b quark mass, Eq. (|TT|). 



6.1 Detector simulation and b-tagging uncertainties 



Biases affecting the jet reconstruction due to the modelling of tracks and clusters were 
estimated as follows. To assess the uncertainty related to the simulation of tracks, the 
resolution of reconstructed track parameters in the Monte Carlo was changed by ±10% [p3fl 
and the analysis repeated. Similarly, to assess the uncertainty related to the simulation 
of the electromagnetic calorimeter, the resolution of this detector was changed by ±10% 
in the Monte Carlo and the analysis repeated. The largest deviation observed was taken 
as the uncertainty due to detector simulation. For all jet finders changing resolution of 
the track parameters by ±10%, i.e. degrading the resolution, gave the largest deviation. 

The effect of the cut on the limited polar angular acceptance of the silicon micro- vertex 
detector was estimated by changing | cos^Thrustl < 0.75 by ±0.05. 

To assess the impact of the b quark selection cut, the analysis was repeated requiring 
for the decay length significance L/ai a minimum value of 8 instead of 5 in the standard 
analysis, which reduced the efficiency to (60.41 ± 0.02) % and lowered the fake tag rate 
to (4.77 ±0.01)%. 

The b tagging efficiency also depends on the mean number of charged particles from 
b hadron decays, their mean lifetime and the branching fraction of Z into bb, which were 
varied within the range given in Tables |5| and || 

The energy spectrum of b and c flavoured hadrons affects both the b tagging and the 
hadronisation. For b quark events the parameter of the Peterson et al. fragmentation 
function f2(J used in the JETSET generator [DJ, et,, was varied from its default value of 
0.0038 by ±0.0010, the range given in [|IR|21],|22|]. Additionally two different fragmentation 



functions were used, Kartvelishvili et al. |32| and Collins and Spiller ||33|| . The largest de- 
viation between the three variations and the standard result was taken as the uncertainty 
due to the modelling of the b fragmentation. For c quark events the parameter of the 
Peterson et al. fragmentation function, e c , was varied from its default value of 0.031 by 
± 0.010 " 



18, 21,22] 
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As shown in Figure [IJ the fraction of b candidates in the data is slightly higher 
than in the Monte Carlo prediction. To estimate the effect of this discrepancy on the 
determination of the b quark mass, the tagging efficiencies e, e 3 (y cut ) were varied by a 
common constant factor. The 1.5% excess of b candidates in data compared to Monte 
Carlo leads to a negligible difference in the mean value for mb, as expected from the small 
uncertainty associated with the variation of i?b- No uncertainty was assigned. 



6.2 Hadronisation uncertainties 



To assess the systematic uncertainties related to the hadronisation process, HERWIG [[35 
was tried as an alternative hadronisation model. It was found that for this generator 
physics involving b quarks is not well described. Among other problems the scaled mean 
energy of weakly decaying b mesons was too low by several standard deviations. Also 
the number of tracks of charged particles found per event was not modelled correctly. 
Therefore no uncertainty was assigned. 

To nevertheless assess the uncertainties related to hadronisation, we altered the main 
parameters of the JETSET model and recalculated the hadronisation correction and the 
resulting value of m^. Beyond the variations affecting the b and c quark fragmentation 
mentioned in the previous section, we altered the value of the b parameter, which affects 
the hardness of the fragmentation function for d,u and s quarks, the width a q of the 
transverse momentum distribution, and the parameter Qq which serves as the cut-off 
for the parton shower. These parameters were varied within their uncertainties quoted 
in Furthermore, the b quark mass inside JETSET was varied by up to ±0.5 GeV in 



0.1 GeV steps. Since a different quark mass significantly affects the details of the hadron 
generation, e.g. emission of soft gluons, which are sensitive to the dead cone effect or 
the formation of b flavoured hadrons, the b quark mass was kept fixed at the JETSET 
default value of 5 GeV throughout the hadronisation process. The variation of the mass 
value was applied only in the calculation of the first gluon radiation probability from the 



b quark which employs the first order matrix element . All these variations are listed 
in Tables E and |6|. 



6.3 Theoretical uncertainties 



Three contributions to the theoretical uncertainty were considered. The coefficients in 
Eq. ( |TTD have uncertainties because of finite statistics used in the Monte Carlo integration 
program. Therefore we altered the coefficients by the uncertainties listed in Table |] and 
reperformed the analysis. This accounts for the uncertainty of the calculation and the 
parametrisation of the mass dependence of the double ratio B 3 . Second, the value of as 
was varied from its world average of 0.1184 by its uncertainty of 0.0031 [51]. Third, the 
renormalisation scale factor in Eq. flip]) was varied by factors of ~ and 2. This last 
variation estimates the impact of neglected higher order terms in the perturbation series. 
The results of these systematic variations are given in Tables || and |6|. 
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JADE 


DURHAM 


JADE E0 


JADE P 


JADE P0 


JADE E 


GENEVA 


CAMBRI. 




result 


0.9752 


0.9532 


1.0657 


0.9938 


1.0140 


1.0649 


0.9817 


0.9717 


statistics 


±0.0048 


±0.0056 


±0.0035 


±0.0040 


±0.0035 


±0.0032 


±0.0094 


±0.0062 


detector simulation 


±0.0047 


±0.0143 


±0.0034 


±0.0076 


±0.0024 


±0.0087 


±0.0218 


±0.0161 


1 cos(0 Thrust ) < 0.75l£$ 


±0.0008 


±0.0011 


±0.0025 


±0.0009 


±0.0013 


±0.0024 


±0.0032 


±0.0024 


L/a L > 8 


±0.0061 


±0.0041 


±0.0038 


±0.0050 


±0.0062 


±0.0089 


±0.0077 


±0.0055 


Ci = 4.955 ± 0.062 


±0.0002 


±0.0010 


±0.0004 


±0.0005 


±0.0003 


±0.0006 


±0.0004 


±0.0006 


r b = (1.564 ±0.014) ps 


±0.0000 


±0.0005 


±0.0003 


±0.0003 


±0.0002 


±0.0003 


±0.0006 


±0.0005 


R h = (0.21751--) 


±0.0015 


±0.0007 


±0.0017 


±0.0013 


±0.0018 


±0.0015 


±0.0004 


±0.0006 


b quark fragmentation 


±0.0088 


±0.0059 


±0.0151 


±0.0114 


±0.0099 


±0.0065 


±0.0102 


±0.0075 


c quark fragmentation 


±0.0007 


±0.0001 


±0.0008 


±0.0003 


±0.0003 


±0.0001 


^0.0015 


^0.0005 


b = (0.52t££)GeV a 


±0.0043 


±0.0013 


±0.0013 


±0.0029 


±0.0037 


±0.0051 


±0.0014 


±0.0027 


* q = (0.40^;^) GeV 


±0.0017 


±0.0012 


±0.0009 


±0.0012 


±0.0011 


±0.0018 


^0.0012 


±0.0009 


Q = (1.90±g : gg) GeV 


±0.0023 


±0.0036 


^0.0010 


^0.0019 


^0.0021 


^0.0010 


±0.0026 


±0.0026 


^b.jETSET = (5.0 ± 0.5) GeV 


±0.0002 


±0.0000 


±0.0011 


^0.0001 


±0.0011 


^0.0014 


±0.0003 


^0.0004 


total systematic uncertainty 


±0.0129 


±0.0166 


±0.0164 


±0.0152 


±0.0129 


±0.0154 


±0.0257 


±0.0192 



Table 5: The systematic uncertainties in the value of the double ratio B 3 . The correlations between the jet finders are taken into 
account. 
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JADE P 


JADE P0 


JADE E 


GENEVA 


mb [GeV 




result 


2.51 


3.20 


2.91 


1.73 


2.58 


2.47 


3.84 


statistics 


±0.42 


±0.24 


to.io 


^0.34 


T-0.12 


T0.07 


±0.55 


detector simulation 


±0.35 


±0.55 


=fo.io 


=F0.72 


=F0.08 
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1 COS^Thrw) < 75 + nn£ 
1 wu^v/ ^ nrust ) — w * 1 0.05 


±0.06 


±0.04 
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T-0.05 
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L/a L > 8 


±0.54 


±0.18 
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decay _ . gr r , q 
"charged — ^- adJ ^ U ' UUZ 
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T , _ (i c;«4 + n (114) ns 
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±0.02 


=fo.oi 
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t-o.oi 


T0.00 


±0.03 


Rh = (0 2175 + nm^ 


±0.11 


±0.03 


^0.05 


t-o.io 


T0.06 


T0.03 


±0.02 


b quark fragmentation 


±0.62 


±0.23 


±0.46 


±1.73 


±0.36 


±0.14 


±0.52 


c quark fragmentation 


±0.05 


^0.00 


^0.02 


T0.02 


t-o.oi 


T0.00 


±0.08 


b = (0.52i££[)GeV 2 


±0.32 


±0.05 
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T-0.13 


tO.h 


±0.07 


cr q = (0.401^3) GeV 


±0.14 
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T0.06 


Qo = (1.90^) GeV 


±0.18 


±0.15 


^0.03 


T-0.15 


T0.07 


T0.02 


±0.14 


^b.jETSET = (5.0 ± 0.5) GeV 


±0.01 


±0.00 


±0.03 


t-o.oi 


±0.04 


T0.03 


±0.01 


total systematic uncertainty 


±0.99 


±0.64 


±0.50 


±1.93 


±0.45 


±0.32 


±1.27 


renormalisation scale 


±0.15 


±0.18 


±0.26 


T0.08 


±0.08 


±0.25 


±0.11 


Aa , Aa l 2 


±0.05 


±0.03 


±0.01 


±0.07 


±0.03 


±0.02 


±0.05 


a s = 0.1184 ±0.0031 


±0.02 


±0.01 


±0.03 


±0.03 


±0.01 


±0.03 


±0.01 


total theoretical uncertainty 


±0.16 


±0.18 


±0.26 


±0.11 


±0.09 


±0.25 


±0.12 


total uncertainty 


±1.09 


±1.20 


±0.56 


±1.86 


±0.52 


±0.41 


±1.50 



Table 6: The uncertainties in the value of the b quark mass at the Z scale. The relative sign between uncertainties for different 
jet finders indicate the sign of their correlation coefficient. 
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0.56 


0.47 
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0.38 
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Table 7: Statistical correlations of the double ratio B 3 between the eight jet finders. The 
coefficients are different from the ones quoted in |K| where 3 and more jet events are used 
to determine the b quark mass. 

7 Combined result 



As can be seen in Table |6| and in Figures || and |6], the individual results of all jet finders 
agree well within their total uncertainty. Hence we now consider a combination of the 
seven determinations. To account for correlations between the seven jet finders, the mean 
was determined using a correlation matrix. This matrix was constructed by dividing both 
data and Monte Carlo into 200 independent subsamplesQ, calculating for each subsample 
and each jet finder the double ratio B3 and determining the correlation between the jet 
finders. This gave the statistical correlation matrix in Table [F[ It differs from that given 
lol as here only 3-jet events are analysed, in contrast to 3 and more jet events in the 
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latter analysis. Using the covariance matrix obtained from the correlation matrix, the 
mean mass value mb of the seven mass values was calculated by minimising 



X 2 = £[5? - B^ the °(m h )] ■ [cov- 1 ]^ ■ [B® - B^ the °(m h )] (12) 



with i = JADE, E, E0, P, P0, DURHAM, GENEVA and cov the covariance matrix of all 
jet finders i. This yielded the result mb(mz) = 2.67 GeV , with a statistical uncertainty 
of ±0.03 GeV and a xVd.o.f. of 100/6. A large xVd.o.f. can be expected since at this 
stage only statistical uncertainties are accounted for. This has been seen also in earlier 
studies fl^ . 



To determine the effect of each of the systematic variations considered in Section ^| on 
the mean value, the diagonal elements of the covariance matrix for the chosen variation 
were added to the statistical covariance matrixQ. With this covariance matrix a new 
mean was calculated with the mass values derived by this particular variation under 
consideration as input values. Its deviation from the standard result was considered as 
the systematic uncertainty due to this variation. The total uncertainty was calculated 
in the same way as for each individual jet finder. Using the same covariance matrix 



3 The inclusive sample contains about one million events and allows for such a fine subdivision 
4 As in |3^| the covariance matrix for a systematic variation of an observable O was constructed by 
assigning the product of two uncertainties SOi ■ dOj from two single measurements i and j to the matrix. 
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Figure 5: For each jet finder the fit to the theoretical prediction, Eq. (|11]), is shown 
together with each measured double ratio B 3 displayed on the y-axis. The b quark mass 
can be read from the x-axis. The band represents the total uncertainty for the double ratio 
and b quark mass. In a) JADE, b) DURHAM, c) JADE EO, d) JADE P, e) JADE PO, f) 
JADE E and in g) GENEVA schemes are shown. No lower bound on the b quark mass 
for the JADE P scheme is shown as this bound reaches the physical limit of a vanishing 
b quark mass. The arrows in a) and c) indicate, how two positively correlated systematic 
uncertainties on B 3 turn into negatively correlated ones in mt,. 



20 



of statistical correlations for all systematic variations assumes that the matrix changes 
only slightly for different variations. The theoretical uncertainty was calculated using the 
weighted mean method. This procedure yields a combined result of 

m h (m z ) = (2.67 ± 0.03 (stat.)^? (syst.) ± 0.19 (theo.)) GeV . (13) 



The uncertainty labelled "syst." includes the detector and the hadronisation terms 
in Tables |5] and |^. The stability of this result was tested by calculating the mean and 
uncertainties using any combination of 2,3,4,5 or 6 out of the 7 jet finders and no large 
deviations from the standard analysis was found. 

To estimate the effect of a non-vanishing c quark mass, the parton level of the light 
quark sample was modified by replacing c quark events by u quark events and repeating the 
analysis. The negligible shift of the mean compared to the standard result was +0.02 GeV, 
which was added in quadrature to the systematic uncertainty. A weighted mean was also 
calculated and gave consistent results. 

Figure |6] shows the b quark mass for the seven individual jet finders together with the 
mean value of the combination and its total uncertainty. 



8 Summary 



The b quark mass was determined at the Z mass scale by comparing the 3-jet rates in b 
and light quark events using seven different jet finders. A deviation of the 3-jet rates in 
tagged b events compared with light quark events was observed. This deviation was used 
to derive the b quark mass by comparing to the theoretical prediction. By minimising 
the x 2 °f th e seven correlated determinations a single result for the b quark mass was 
determined to be 

m b (m z ) = (2.67 ± 0.03 (stat.)^? (syst.) ± 0.19 (theo.)) GeV . (14) 

Our final result is shown in Figure |7|. Also shown in this figure is the average of the b quark 
mass at the scale of the b quark mass itself, mf DG '(mb) = (4.2 ± 0.2) GeV. as compiled 
in [[3] . This average has been derived from measurements at the bb production threshold 
and from b hadron masses. Also shown are other measurements of m^mz) by DELPHI || , 



ALEPH H and Brandenburg et al. using SLD data ||10|| . These determinations at the 
Z mass scale yielded mass values in the range of 2.67 GeV to 3.27 GeV. All results are 
in good agreement with each other. The total uncertainty on the b quark mass in this 
analysis is smaller than for previous measurements at the Z mass scale, see Figure [f[ The 
solid curve in Figure ^ shows the QCD expectation of a scale dependent b quark mass in 
the MS renormalisation scheme, using the world average value of ag(m|) = 0.1184±0.0031 



[5TJ. Evolving our result of mb(mz) down to the b quark mass scale itself gives 



m h (fn h ) = (3.95 ±0.04 (stat.)^ (syst.) ±0.28 (theo.)) GeV 
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Figure 6: The result of the combination of the seven individual measurements is shown 
by the vertical line, the hatched band displays its total uncertainty Also shown are the 
mass values with their uncertainties for each individual jet finder. The inner error bars are 
the statistical uncertainty. The arrows depict the statistical and systematic uncertainties 
added in quadrature, and the outer error bars depict the total uncertainty. 
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Figure 7: Results for rrlb(mz) from this analysis and at the scale % from the compilation 
in [[| , together with other determinations at the Z mass scale from DELPHI || , ALEPH |J 
and Brandenburg et al. . The solid curve shows the theory expectation for the running 
of the b quark mass, Eq.(@), using a s (m 2 z ) = 0.1184 and m b PDG) (m b ) = 4.2 GeV §. 
The dotted lines display the total uncertainty on the b quark mass run from production 
threshold to the Z pole. This total uncertainty includes the uncertainty on the mass 
itself, m b PDG) (m b ) = (4.2 ± 0.2) GeV and the uncertainty on a s (m|) = 0.1184 ± 0.0031. 
The error bars of the data points show the statistical, the systematic and the theoretical 
uncertainties added in quadrature. For displaying purposes the four measurements at the 
Z pole have been separated. 
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which is in agreement with (4.2 ± 0.2) GeV j3|. 

We have also compared our result for the b quark mass at the Z mass scale with the 
combined value at production threshold, yielding 

mf DG) (m b ) - m b (m z ) = (1.53 ± 0.39) GeV , (15) 

which is different from zero by 3.9 standard deviations, confirming the running of the b 
quark mass in the MS renormalisation scheme as predicted by QCD. 
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